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Abstract 

A feasibility study on the use of shape memory 
alloys (SMA) for suppression of the random response 
of composite panels due to acoustic loads at elevated 
temperatures is presented. The constitutive relations 
for a composite lamina with embedded SMA Rbers are 
developed. The finite element governing equations and 
the solution procedures for a composite plate subjected 
to combined acoustic and thermal loads are presented. 
Solutions include: 1) critical buckling temperature, 

2) flat panel random response, 3) thermal postbuck- 
ling deflection, and 4) random response of a thermally 
buckled panel. The preliminary results demonstrate 
that the SMA fibers can completely eliminate the ther- 
mal postbuckling deflection and significantly reduce 
the random response at elevated temperatures. 

Introduction 

Future high-performance aircraft, such as the high 
speed civil transport (HSCT), will subject the external 
skin panels to significant acoustic and thermal loads. 
The new vectored thrust propulsion systems on the 
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YF-22 Advanced Tactical Fighter ATE provide short 
take-off and landing capability and increased maneu- 
verability. However, the vectored thrust also subjects 
portions of the aircraft structure to intensified acous- 
tic and thermal loads. Acoustic fatigue thus emerges 
as one of the critical considerations in the design and 
development of high speed flight vehicles k 

A significant amount of research in the last few 
years has greatly advanced the held of control of flex- 
ible structures through the use of smart materials^. 
The held of adaptive or intelligent structures involves 
a broad range of technologies that enable the realiza- 
tion of structural systems that are able to sense and 
to control their own behavior, such that the range 
of operational performance for these materials may 
be extended over conventional materials and struc- 
tures. Smart materials have the ability to change 
shape, stiffness, natural frequency, damping, and other 
mechanical characteristics in response to changes in 
temperature, electric held, or magnetic held. The 
most common smart materials are shape memory al- 
loys, piezoelectric materials, magnetostrictive materi- 
als, electrorheological fluids and magnetorheological 
fluids. 

Shape Memory Alloy Hybrid Composites 

Shape memory alloys exhibit a characteristic 
phase transformation from martensite to austenite, 
initiating at the austenite start temperature © and 
asymptotically ending at the austentite finish tempera- 
ture Tj . A shape memory alloy in the low temperature 
martensitic condition (T < ©), when plastically de- 
formed and the external stresses removed, will regain 
its original (memory) shape when heated. For exam- 
ple, strains of typically six to eight percent can be com- 
pletely recovered by heating the nickel-titanium alloys 
(Nitinol) above the austenite finish temperature. The 
transformation temperatures can be altered by chang- 
ing the composition of the alloy. In addition, when 
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Nitinol is heated, the Young’s modulus increases three 
to four times and the yielding strength also increases 
approximately ten times^^. 

The materials, referred to in this study as SMA 
hber-reinforced hybrid composites, are conventional 
advanced composite materials (such as graphite-epoxy) 
that contain embedded SMA Rbers having the same di- 
rection as the graphite fibers. The memory effect of 
the SMA fibers may be put to use by applying an ini- 
tial elongation to them during fabrication of the lami- 
nate. Thus, once the laminate is heated above Tg, the 
SMA fibers will try to recover the initial strain and in- 
duce tensile inplane forces in the laminate. Since the 
expected aerodynamic heating on the skin panels of 
a HSCT during supersonic cruise will generate panel 
temperatures well above the austenite finish tempera- 
ture, SMA will be ideal for such an application. The 
overall stiffness of the hybrid composite panel will be 
increased due to: 1) the increase of the Young’s modu- 
lus of the SMA fibers by a factor of three or four, and 
2) the internal tensile inplane forces induced in the 
panel from the recovery of initial strains of the SMA 
fibers. Therefore, the root-mean-square (RMS) maxi- 
mum deflection and RMS maximum strain, due to ran- 
dom pressure loading, will be reduced at elevated tem- 
peratures as compared to composite panels without 
embedded SMA fibers. This has been demonstrated in 
a preliminary investigation on passive control of ran- 
dom response of SMA hybrid composite plates using 
the classic analytical continuum method^. 

Shape memory alloys have been applied as ac- 
tuators for active control of buckling of beams* and 
shape control of beams®. R is also being studied for 
use in active vibration control of beams^®di and large 
space structures^®. Active vibration control of flexible 
linkage mechanisms using SMA hber-reinforced com- 
posites has been investigated by Venkatesh et al.^*. 
Acoustic transmission and radiation control through 
the use of SMA fibers in a hybrid composite was 
presented by Liang and Anders^"^d5 all of these 
investigations*^^*, the SMA fibers or strips were heated 
by applying an electrical current with some control de- 
vice. Since it was assumed that the composite would 
not be appreciably heated, thermal expansion effects 
were not included. However, the heat source for the 
present study is assumed to be due to aerodynamic 
heating. Therefore no control device is needed, but 
thermal expansion effects are important since the en- 
tire laminate is heated and thermal buckling and post- 
buckling considerations may be necessary. 


the graphite fibers. Figure 1. The 1 and 2 directions, 
or principal material directions, are parallel and per- 
pendicular to the fiber direction respectively. The con- 
stitutive relations for such a lamina can be derived fol- 
lowing the engineering approach found in many com- 
posite materials texts. A derivation may be found in 
the appendix and the resulting constitutive relations, 
in principal material coordinates, are given by Eqs. 
(A17) and (A18). For a general k-th layer with an ori- 
entation angle 9, the stress-strain relations, Eqs. (A17) 
and (A18), become 



( 2 ) 


The resultant force and moment vectors of the SMA 
hber-reinforced hybrid composite plate are defined as 


{N,M} 



or 
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where the laminate stiffness [A*], [B*] and [D*] 

and the recovery inplane force and moment vectors 
{#*} and {M*} are all temperature dependent. The 
inplane strain and curvature vectors are defined from 
the von Karman strain-displacement relations 
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Finite Element Formnlation 


Consider a thin composite lamina, for example 
graphite-epoxy, having arbitrary orientation angle 6 
and embedded SMA fibers in the same direction as 


A limited number of investigations on the struc- 
tural response of panels subjected to combined acous- 
tic and thermal loads exist in the literature. Seide and 
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Adami^® were the first who studied large deflection 
random response of a thermally buckled beam. The 
well-known classic Woinowsky-Krieger large amplitude 
beam vibration equation was used. Galerkin’s method 
and time domain numerical simulation were then ap- 
plied to obtain the random response. More recently, 
the Galerkin/numerical simulation approach was ap- 
plied to simply supported metal and orthotropic com- 
posite rectangular plates by Vaicaitis et al.^^d*^ 'j'jjg 
classic von Karman large-deflection plate equations in- 
cluding uniform temperature and orthotropic property 
effects were employed. Lee^® extended that work to in- 
clude nonuniform temperature distributions. He stud- 
ied isotropic rectangular plates with either simply sup- 
ported or clamped edges. The Galerkin/equivalent lin- 
earization method^® was used. The classic continuum 
approaches thus have been limited to simple beams^® 
and isotropic or orthotropic rectangular plates^^^^®. 

For over three decades, the Rnite element method 
has been the predominant approach for structural me- 
chanics of complex structural geometry. However, 
there are few studies where the random response of 
structures subjected to combined acoustic and ther- 
mal loads is involved. Locke and Mei^^’^^ extended the 
finite element method for the first time to beam and 
plate structures under combined thermal and acous- 
tic loads. The thermal load considered was a steady- 
state temperature distribution AT(x,y). Ghen and 
Mei^3 recently refined that finite element formulation 
for nonlinear random response of structures by con- 
sidering the acoustic pressure and thermal load to be 
applied simultaneously. 

The governing equations of motion can be derived 
for a SMA fiber-reinforced hybrid composite panel sub- 
jected to combined thermal and acoustic loads through 
the use of a variational principle. These equations may 
be written in the form 
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where [M] and [A'] are the system mass and stiff- 
ness matrices; [A'*] and [Knat] are the geometric 
stiffness matrices due to the recovery stress a* (or 
{A*}) and thermal inplane force vector {Aat}, re- 
spectively; [Al] and [A2] are the first and second or- 
der nonlinear stiffness matrices which depend linearly 
and quadratically upon displacement {Lh}, respec- 
tively; {A(t)| is the random acoustic excitation load 
vector, {Pat} is the thermal load vector, and 
is the SMA recovery force vector. The subscripts h 
and m denote bending and membrane components, 
respectively; and the subscripts B, Nm and NB 
indicate that the corresponding stiffness matrix is 
due to the laminate stiffness [B], membrane forces 
{A„}(= [A][e°^}) and {As}(= [B]{k}), respectively. 
The stiffness matrices [A'], [Al] and [A2] and the re- 
covery force vector are also temperature depen- 

dent. Detailed derivations of the equations of motion 
and expressions for the element matrices and load vec- 
tors are referred to reference 21. The solution proce- 
dures developed in references 21-23 are to be extended 
to obtain the critical buckling temperature, thermal 
postbuckling deflection, and linear RMS maximum de- 
flection of SMA hybrid composite panels. 


Solution Procedures 

Four different types of analyses are required for 
sonic fatigue design of SMA hybrid composite panels. 
They are: (1) thermal buckling analysis, (2) random 
vibration analysis, (3) thermal postbuckling analysis, 
and (4) random vibration analysis of thermally buck- 
led panels. Each analysis is discussed briefly in the 
following. 

Thermal Buckling Analysis 

As the first step, a thermal buckling analysis is 
carried out to determine the critical buckling temper- 
ature Agr = ATgr + To , where T„ is the reference or 
ambient temperature. Equations governing the ther- 
mal buckling are obtained by neglecting the inertia, 
nonlinear, bending-membrane coupling, and bending 
force terms from the incremental form of Eq. (5) as 

([Ai] + [a:] - [Knat] + [AlAr™]){Am} = 0 (7) 

and 

[K^]{W^} = {P^at} - {A^,} (8) 

The first-order nonlinear stiffness matrix [Aljvm] is 
linearly dependent upon the inplane displacement 
{WGj} which can be expressed at a given AT as 

= [A™]-'({A„at} - {A^,}) 

= {Wm}^-{Wm}„ 
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where {Wm}i = [Km] ^{PmAr} is linearly dependent 
upon AT and {Wm}g = [Km]~^ {Pma} i® ^ constant 
vector. Thus, the matrix [Atljvm] is the sum of two 
matrices: the hrst matrix [Atljvm]i is evaluated with 
{Wm}i and the second matrix [Atljvm]o is evaluated 
with {Wm}„ as 

[iVliv™] = [tVliv™]i - [tVliv™]„ (10) 

The matrices [A'jvat] and are linearly de- 

pendent upon AT. Therefore, Eq. (7) can be ex- 
pressed as an eigen-problem 

{[Kb] + [K*] - [NlNm]o){<f>] = K[Knat] 

The critical temperature change corresponds to 
the lowest value of A and is given by AT^^ = 
AAT, and {(f)} is the corresponding buckling mode. 
Since the matrices [K^], [K*], [Atljvmlo and 
in Eq. (11) are all temperature dependent, an itera- 
tive scheme is implemented to achieve convergence for 
AT„ (A = 1.000 ± .001). 

Random Vibration Analysis 

Once the critical buckling temperature is deter- 
mined, the next step will depend on the panel tem- 
perature for the particular application of concern. Eor 
panel temperatures greater than the critical tempera- 
ture (T > Ter), a thermal postbuckling analysis has to 
be performed, followed by a random vibration analysis 
of the thermally buckled structure. 

Eor panel temperatures less than the critical tem- 
perature (T < Ter), the panel is flat. A random vibra- 
tion analysis including the thermal and SMA effects 
is performed. The response results include frequen- 
cies and modes of vibration, RMS maximum deflec- 
tion, and power spectral density (PSD) functions of 
maximum deflection. The equations of motion for a 
SMA hybrid composite panel subjected to a combined 
acoustic excitation and thermal load can be obtained 
from Eq. (6) as 

[M]{w] + [K]{W] = {Pit)} (12) 

where 

[K] = [K] + [I<:] - [Knat] + [tvi] (13) 

where the only nonzero entries in [Nl] are due to 
[NInjjj] resulting from the inplane deflection, Eq. (9). 

The natural frequencies ujr and mode shapes {(fir} 
of vibration are obtained from the eigenvalue problem 

UJ^r[M][(}r} = [mc)r} (14) 


A set of uncoupled modal equations with reduced 
degrees-of-freedom can thus be obtained from Eq. (12) 
as 


with the truncated modal transformation 
N 


{IE} = '^{<j)r}qr = [(f]{q} 


(16) 


In Eq. (15), a modal damping term has been included 
where C is the modal damping ratio and the modal 
mass and force are 


rUr = {()>r} [Mb]{(j>r} 

fr = {(t)r}^{Ph(t)}/mr 


(17) 


The random modal response can be easily determined 
from Eq. (15) for the case of a spatially uniform white 
noise random load p(t) with spectral density Sp or Gp 
as 


/ OO 
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with Gpif) = iirSpiuj) for w = 2irf (/ > 0) and 
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The cross-correlation of modal amplitudes from Eq. 
(18) is 


E[qrqs] = 


Gp(f)(Cr‘^r + Cs^s) 


nirmeiioJ^ -w'j) +A[CrCsUlr^s{i^? + ^I) 


+ (Cr + Cl)‘^l‘^l]] 


(20) 

and the correlation matrix of the system nodal ampli- 
tude becomes 


E\{W}{Wy =m {q]{q} 


(21) 


Thermal Postbuckling Analysis 

Eor panel temperatures greater than the critical 
temperature (T > Ter), postbuckling analysis is per- 
formed to determine the thermal deflection. The gov- 
erning equations for thermal postbuckling can be ob- 
tained from Eq. (6) as 

([4t-1 + [K{] — [A'jvat] + 

= {Aat}-{a;i 

( 22 ) 
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One effective approach involves using the Newton- 
Raphson iterative method for solving time-independent 
nonlinear problems. Thus, for the j-th iteration the in- 
cremental form of Eq. (22) can be written as 

[K,,n],{AW}^^,={AP}^ (23) 

Then [A"tan]j_|_i and are updated by using 

{1Tat}j_|_i = {Wat}^ + {AW]-j^-^ (24) 

The solution process seeks to reduce the imbalance 
load vector {AA}, and consequently {AVE}, to a spec- 
ihed small quantity. The tangent stiffness matrix and 
the imbalance load vector are 

[A'tanlj = [K] + [A'*] - [A'jvat] + + [^^2at]j 

(25) 

and 

{AA}^. = {Pat} ~ {a;| - (^[A] + [A':] - [AAat] 

+ ^[NlAT]j + ■^[N‘2AT]j^{WAT}j 

(26) 

where the subscript AT denotes that the nonlinear 
stiffness matrix is evaluated using {VEat}- The linear 
buckling mode shape, Eqs. (9) and (11), normalized 
by a scale factor, is usually taken to be the initial trial 
solution for Eq. (23). 

Random Vibration Analysis of Bnckled Panels 

The random vibration analysis of a thermally 
buckled panel is performed to determine the frequen- 
cies and modes of vibration about the buckled po- 
sition, the RMS maximum deflection, and the PSD 
functions. The equations of motion for a thermally 
buckled SMA hybrid composite panel can be obtained 
from Eq. (6) and ref. 24 as 

[M]|it| -1- ([A] -I- [K*] - [Knat] + [AIat] 

+ [N2at]){W} = {P{t)} 

where the subscript AT denotes that the correspond- 
ing stiffness matrix is evaluated using the thermally 
buckled deflection {VEat}- The influences of re- 
covery stress a*, thermal force {Vat} and thermal 
postbuckling deflection are given by the presence of 
[A*], [AAat], [AIat] and [A2 at] in Eq. (27). An 
examination of Eqs. (25) and (27) reveals a novel ap- 
proach of the present Rnite element formulation. That 
is, the stiffness matrices in Eq. (27) need not to be 
computed and assembled. The sum of those stiffness 


matrices is exactly the converged tangent stiffness ma- 
trix in thermal postbuckling analysis. Thus, Eq. (27) 
is simply 

[M]{EE} + [Atan]{fff} = {A(t)} (28) 

The natural frequencies and mode shapes of vi- 
bration about the thermally buckled position {VEat} 
are determined from 

= [AAn]{<(>J (29) 

Equation (27) thus becomes a set of uncoupled modal 
equations 

Ss + 2Cs^sSs + Ss = /s ) s = l,2,...,A (30) 

with the truncated modal transformation 

N 

{EE} = ^ {<l)s}qs = [<!>]{q} (31) 

S = 1 

The determination of the random response follows 
from Eqs. (17)-(21). 

Results and Discussion 

The results shown in this paper were generated for 
a SMA hybrid composite laminate, where the graphite- 
epoxy composite was treated as the matrix. The fol- 
lowing material properties were used in the analysis: 

Graphite-Epoxy 


El 155 GPa (22.5 Msi) 

Ea 8.07(1.17) 

Gi2 4.55 (0.66) 

Vi2 0.22 

P 1550.07 Kg/m^ (0.1458x10“^ Ib-s^/in."^) 

-0.07x10-® /°C (-0.04xl0-V°E) 

«2 30.1x10-® (16.7x10®) 

SMA-Nitinol 

T, 37.78 °C (100 °E) 

Tf 62.78 (145) 

E Erom Eig. A3 
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G 24.86 (3.604) for T < 

25.6 (3.712) for T > Tf 

V 0.3 

P 6450 Kg/m^ (0.6067x10“^ Ib-s^/in."^) 

a 10.26x10-® /°C (5. 7 x10-V°F) 

Table 1 presents the results for a simply supported 
rectangular laminate with a graphite-epoxy matrix. 
The dimensions of the laminate are 15x12x0.048 in. 
and the layup is symmetric (0/45/-45/90)g. Results 
are presented for a laminate with no SMA Rbers and 
for three SMA initial strain values (3%, 4% and 5%) 
for each nonzero SMA volume fraction. The head- 
ings indicate SMA volume fraction, SMA fiber initial 
strain, two critical temperature changes, normalized 
maximum thermal deflection, and normalized maxi- 
mum RMS deflection. 

Note that when SMA fibers are present in the lam- 
inate, there are at least two critical buckling tempera- 
tures (ambient temperature 70°F). One corresponds to 
temperatures less than the austenite start temperature 
Tg , and the other corresponds to temperatures greater 
than the austenite finish temperature Tj. It can be 
seen that the addition of SMA fibers slightly reduces 
the critical buckling temperature below Tg . This is 
caused by the relatively low modulus of the SMA fibers 
in relation to the graphite epoxy matrix they are re- 
placing. However, with the addition of only 10% SMA 
fibers at 3% initial strain, the critical buckling temper- 
ature above Tj is increased drastically. Some of the 
reported critical buckling temperatures are unrealistic 
for a graphite-epoxy matrix since epoxy composites 
are generally considered to have a maximum service 
temperature of around 300° F. However, these results 
are qualitatively indicative of the shape memory effect 
and allude to the potential use of SMA fibers in high 
temperature matrices such as metals, carbon-carbon, 
etc. 

It can be seen that the thermal deflection W^r^/h 
is attenuated in all cases and even eliminated in most. 
The themal deflection is eliminated for the cases with 
10% SMA fibers at 5% initial strain and all greater 
volume fractions since the critical buckling tempera- 
ture has been shifted well above the panel tempera- 
ture (300°F). This effect is beneficial due to the fact 
that the optimal structural configuration may be main- 
tained and the possibility of snap-through motions is 
eliminated. 

Although the dynamic response Wmax/h of panels 
with some of the lower SMA volume fractions and ini- 
tial strains is actually greater than that for the panel 
without SMA fibers, significant dynamic response at- 


tenuation may be achieved with SMA volume frac- 
tions of thirty percent or more. The increase in the 
dynamic response, observed for some of the cases, is 
due to the relatively low modulus of the Nitinol alloy, 
even at high temperatures and another factor which 
will be discussed subsequently. Finally, the case of 
10% SMA and 4% initial strain exemplifies the im- 
portance of choosing the SMA parameters carefully, 
as indicated by the excessively large RMS maximum 
deflection. This particular choice results in a critical 
buckling temperature of 295.06°F, which is within 2% 
of the panel temperature (300°F). Thus, the panel is 
very compliant since it is close to buckling. 

The variation of the RMS (Wmax/h) with sound 
spectrum level (SSL) is shown in Figure 2 for four 
cases: no SMA fibers at 70°F, 10% SMA fibers and 
5% initial strain at 300°F, no SMA at 300°F, and 30% 
SMA and 5% initial strain at 300°F. Note that the 
panel without SMA fibers at 300°F exhibits signifi- 
cantly reduced dynamic response because of the addi- 
tional stiffness due to the fact that the panel is ther- 
mally buckled. Adding 10% SMA fibers to the 300°F 
panel increases the dynamic response since the recov- 
ery forces induced are not sufficient to overcome the 
modulus deficiency of the SMA fibers and the loss of 
stiffness due to the fact that the panel is no longer ther- 
mally buckled. Finally, addition of 30% SMA fibers 
provides ample recovery forces to significantly reduce 
the dynamic response, even more than the thermally 
buckled case. 

The variation of the RMS (Wmax/h) with temper- 
ature for a laminate with no SMA fibers in comparison 
to one with 30% SMA and 5% initial strain is shown 
in Figure 3. Note that they both exhibit a peak, albeit 
beyond the linear range (Wmax/h<l), at their respec- 
tive low critical temperatures due to extremely low 
stiffness at this point of instability. In addition, the 
response of the laminate with SMA fibers has a peak 
at 100°F since the SMA modulus has a local minimum 
at that temperature. 

The power spectral density of the maximum de- 
flection at 100 dB is shown in Figures 4 — 6 for three 
cases: no SMA at 70°F; 30% SMA and 5% initial strain 
at 70°F; and 30% SMA and 5% initial strain at 300°F, 
respectively. At low temperature, the fundamental fre- 
quency of the panel with SMA fibers is lower than that 
for the panel with no SMA since the SMA reinforced 
panel is less stiff at low temperatures and more mas- 
sive. At 300°F, the modes of the SMA Rber-reinforced 
panel exhibit significant amplitude reduction and have 
shifted to much higher frequencies. 
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Conclusions 

The finite element method has been successfully 
implemented to analyze the thermal buckling, ther- 
mal postbuckling, and random response of SMA Rber- 
reinforced hybrid composite panels. The constitu- 
tive relations, including the thermal and shape re- 
covery effects, have been derived and used in the fi- 
nite element formulation to model the mechanics of a 
SM A/graphite-epoxy panel. 

It was found that panels with SMA fiber reinforce- 
ment have at least two critical buckling temperatures; 
one below the austenite start temperature Tg , and one 
above the austenite finish temperature Tj. The ex- 
tremely high critical buckling temperatures predicted 
for SMA volume fractions greater than 10% indicate 
the possible suitability of SMA fibers in high temera- 
ture composites. 

The thermal postbuckling deflection can be to- 
tally eliminated, for panel temperatures greater than 
Tf, with as little as 10% SMA fibers. The property 
could be useful in practical applications by maintain- 
ing optimal aerodynamic configuration and eliminat- 
ing snap-through motions. 

The RMS maximum deflection actually increases 
for some SMA volume fractions because of two fac- 
tors: the modulus of the SMA fibers is relatively low, 
even at high temperatures; and the recovery forces in- 
duced by small volume fractions of SMA fibers is sufh- 
cient to prevent thermal buckling (for Tf<T<Tcr) but 
insufhcient to overcome the loss of the buckling stiff- 
ness. However, significant reductions in the dynamic 
response are realizable with SMA volume fractions of 
30% or more. 
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Figure 1: SMA Rber-reiforced hybrid compostite lam- 
ina. 
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Table 1: Results for simply supported (0/45/-45/90)s 
panel (ambient temperature 70°F). 



Figure 2: RMS maximum deflection versus sound 
spectrum level. 
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Figure 3: RMS maximum deflection versus tempera- 
ture. 
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Figure 4: PSD of Wmax for panel without SMA at 
70°F. 
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Figure 5: PSD of Wmax for 30% SMA, 5% initial strain 
at 70°F. 
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Figure 6: PSD of Wmax for 30% SMA, 5% initial strain 
at 300°F. 

Appendix 

Hybrid Lamina Constitntive Relations. 

A representative volume element of a SMA fiber- 
reinforced hybrid composite lamina is shown in Figure 
7. The element is taken to be in the plane of the 
plate. The composite matrix, for example graphite- 
epoxy, has the principal material directions 1 and 2 
and the SMA Rber is embedded in the 1-direction. 

In order to derive the constitutive relation for the 
1-direction, it is assumed that a stress <7i acts alone 
on the element (<72 = 0) and that the SMA fiber and 
composite matrix are strained by the same amount, 
€i (i.e., plane sections remain plane). The 1-direction 
stress-strain relation of the SMA fiber can be described 
as 

au = K^i+^: , T>T, (A1) 

or 

,71, =A,*(ei-a,AT) , T < T, (A2) 

where T, is the austenite start temperature and a, is 
the thermal expansion coeffecient. The Young’s mod- 
ulus E* and the recovery stress a* are temperature de- 
pendent, indicated by the superscript (*). The recov- 
ery stress a* is also dependent on the initial strain 
For Nitinol, a* and E* can be determined from Figures 
8 and 9, respectively"^. Similarly, the one-dimensional 
stress-strain relation in the 1-direction for the com- 
posite matrix can be expressed as 

^ImAT) (-^^) 

The resultant force in the 1-direction (ct 2 = 0) is 
distributed over the SMA fiber and composite matrix 
and can be written as 

(A4) 
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where {(Ti,Ai), {(Tis,As), and (crim,Aim) are the 
(stress, cross sectional area) of the entire element, 
SMA hber, and composite matrix, respectively. Thus, 
the average stress ai is 


ai = ai,Vs + (TimVm (AS) 

where Vg = Ag/A\ and v^n = AmjAx are the volume 
fractions of SMA and composite matrix, respectively. 
When T > Tg, the SMA effect is activated and the one- 
dimensional stress-strain relation in the 1-direction 
becomes 


(Ti = {E*g€i + (T*)vg + Eim{ei - aim/\T)v,n 

= EleiA a*Vg - Eimai^Vm^T 


(A6) 


where 

E*i = EimVm + E*gVg (A7) 

When T < Ag, the SMA effect is not activated and 
the stress ai is 


(Ti = - {E*agVg + Ai™ai™t;™)AT 

= A((ei - aiAT) 

where 

Qfi — 

A similar constitutive relation may be derived for 
the 2-direction by assuming that the applied stress 
(72 acts upon both the hber and the matrix (ai = 0). 
Thus, the one-dimensional stress-strain relations in 
the 2-direction for the SMA hber and the composite 
matrix become 


(A8) 

(A9) 


a2g = a2 = E*g{e2g-agAT) (AlO) 

and 

<^2m = = E2m{<^2m ~ (AH) 

respectively. The recovery stress does not appear in 
Eq. (AlO), since the SMA hber initial strain and 
recovery stress a* are considered to be a 1-direction 
effect only. 

The total elongation is due to strain in the com- 
posite matrix and the SMA hber and may be written 
in the form 


Ai€ 2 — Ami2m -\- Ag€2g (A12) 


Thus, the total strain becomes 


Since (T 2 = (^2 — a 2 AT), Eqs. (AlO) and (All) 

may be substituted into Eq. (A13) to give 


£2 — + (X2AT 

Et 


^2 

(^2Vs 

E* 


0'2Vn 

E2.rr 


A {agVg A a2m'Vm)^T 


(A14) 


Therefore, the modulus and thermal expansion coefh- 
cient in the 2-direction become 


and 


^ E2mEl 


®2 — 0!2m'Vm + 


(A15) 

(A16) 


Expressions for the hybrid composite Poisson’s ratios 
and shear moduli follow from similar derivations. 

The constitutive relations for a thin composite 
lamina with embedded SMA Rbers can be derived 
using a similar engineering approach to give 


{ <Tl ) 

'Qh 

0)2 

0 

1^^ ] 

.2 = 

Qi2 

022 

0 

1 "" 

1 E2 I 

0 

0 

0(36 . 

1 712 J 



AT 


(A17) 


T>Tg 


and 



= [01 1^1 eS jaj ^AT| , T < Tg 

(A18) 

where [Q]^ and [Q*] are the reduced stiffness matrices 
of the composite matrix and the composite lamina, re- 
spectively. The [Q*] matrix is temperature dependent 
and is evaluated using the previously derived relations 
as 


(E*i,12i 2) = {Eim,Vl2m)Vm +{E*g,Vg)Vg (A19) 


and 


{E*2,G\2) = 


{E2mE* , Gi2mG*) 


[{E2m, Gi2m)Vg A (E* , G*)Vm\ 


(A20) 


£2 — £2m'l’m + i2sV, 


(A13) 


10 


where the As are Poisson’s ratios and the G’s are 
the shear moduli. The thermal expansion coefhcients 
ai and o ;2 are derived in Eqs. (A9) and (A16). 
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Figure 7: Representative volume element for SMA 
fiber-reinforced hybrid composite lamina. 
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Figure 8: SMA recovery stress versus temperature and 
initial strain, ref. 4. 
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Figure 9: SMA modulus versus temperature, ref. 4. 
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